Near-Infrared Survey of the GOODS-North Field: Search for 
Luminous Galaxy Candidates at z > 6.5 *t 

Nimish P. Hathi 1 ' 2 , Bahrain Mobasher 2 , Peter Capak 3 , Wei-Hao Wang 4 , 

Henry C. Ferguson 5 

nhathiOobs . carnegiescience . edu 
ABSTRACT 

We present near-infrared (NIR; J & K s ) survey of the Great Observatories 
Origins Deep Survey-North (GOODS-N) field. The publicly available imaging 
data were obtained using the MOIRCS instrument on the 8.2m Subaru and the 
WIRCam instrument on the 3.6m Canada- France-Hawaii Telescope (CFHT). 
These observations fulfill a serious wavelength gap in the GOODS-N data - 
i.e., lack of deep NIR observations. We combine the Subaru/MOIRCS and 
CFHT/WIRCam archival data to generate deep J and K s band images, covering 
the full GOODS-N field (~169 arcmin 2 ) to an AB magnitude limit of ~25 mag 
(3<t). We applied z 850 -band dropout color selection criteria, using the NIR data 
generated here. We have identified two possible Lyman Break Galaxy (LBG) 
candidates at z > 6.5 with J < 24.5. The first candidate is a likely LBG at 
z ~ 6.5 based on a weak spectral feature tentatively identified as Lya line in the 
deep Keck/DEIMOS spectrum, while the second candidate is a possible LBG 
at z ~ 7 based on its photometric redshift. These ^so-dropout objects, if con- 
firmed, are among the brightest such candidates found so far. At z > 6.5, their 
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star formation rate is estimated as 100-200 M yr" 1 . If they continue to form 
stars at this rate, they assemble a stellar mass of ~5xl0 10 M after about 400 
million years, becoming the progenitors of massive galaxies observed at z ~ 5. 
We study the implication of the z 850 -band dropout candidates discovered here, 
in constraining the bright-end of the luminosity function and understanding the 
nature of high redshift galaxies. 

Subject headings: catalogs — galaxies:formation — galaxies: evolution — galaxies: 
high-redshift 



1. Introduction 



In recent times, the high redshift frontier has started to probe the epoch of reioniza- 
tion because of extensive observations at Near-Infrared (NIR) wavelengths. This extension 
of multi-wavelength galaxy surveys to obtain deep images at observed NIR has provided a 
new prospect for the study of galaxies at high [z > 6.5) redshifts. Such observations are 
vital to identify star-forming high redshift objects by sampling their rest-frame UV light, 
and for estimating their star formation properties. On the other hand, the combination 
of Spitzer /IRAC and NIR observation s can be used to id entify massive, evolved Balmer 



Break Galaxies (BBGs) at z > 5 (e.g.. IWiklind et al.l 120081 ). Furthermore, relative insensi- 



tivity at these wavelengths to the stellar population mix, dust and K-correction, specially at 
intermediate (z ~ 2) redshifts is an added bonus. 

Among the most extensive datasets avai lable in extragalactic a stronomy, is the Great 
Observatories Origins Deep Survey (GOODS: iGiavalisco et al.ll2004ai ). This consists of imag- 
ing of two fields (GOODS-North and GOODS-South), by the Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) in four optical bands (-B435, Vqoq, 2775, ^85o)- Given 
the need for NIR data in any study of evolution of galaxies, the GOODS-South (GOODS-S) 
field h as been observed ext ensively at NIR wavelengths, using the 8.2m Very Large Telescope 
(VLT; iRenzini et all l2003h in the J, H and K s bands, but the GOODS-North (GOODS- 
N) field has not had such full-field deep coverage at these wavelen gths. In the last few 
years, the Wide-field InfraRed Camera (WIRCam; iPuget et al.l 120041 ) on the 3.6m Canada- 
France- Hawah^Tj^ks^ojDe^ Infrared Camera and Spectrograph 
(MOIRCS; llchikawa et al.l l2006t ISuzuki et al.l 120081 ) on the 8.2m Subaru telescope have in- 
vested large amounts of observing time to obtain NIR observ ations of the GOODS-N field. 
The WIRCam K s data wer e published by I Wang etal.l hold ) and the MOIRCS J, H, and 
K s data were published by iKajisawa et al.l (120111 ). In both works, the images and catalogs 
were made publicly available. Various studies have used subsets of these observations for 
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specific science goals. For example, iKajisawa et al . (2006) used MOIRCS data to investigate 
the number counts of Distant Red Galaxies, while llchikawa et al.l (120071 ) used MOIRCS K s - 
band selecte d galax ies to measure clustering properties of low-mass galaxies at 1 < z < 4. 
Bundy et al.l (120091 ) used NIR observations from both GOODS fields to investigate the de- 
pend ence on stellar mass an d galaxy type of the close pair fractions and implied merger rate, 
while iBouwens et al.l ( 120081 ) used the deepest NIR observations from both fields to search for 
star-forming galaxies at z > 6.5 and constrain their rest-frame ultraviolet (UV) luminosity 
functions. These studies show the significance of NIR observations in such well-studied fields. 
The full coverage of the GOODS-N region is essential not only to expand these studies, but 
also to accomplish larger number of additional science goals. Despite such studies, there is 
limited information in literature of the observational details and sensitivities of the NIR data 
for GOODS-N field, and details about source detection. 

The high sensitivity of the HST /Wide Field Camera 3 (WFC3) has made it possible to 



identif y faint Lyman break galaxies (LBGs) at z > 6.5 (e.g..lOesch et al.ll2010t IBouwens et al. 



2010aU bl; iFinkelstein et al.l I201CH : lYan et al.l |2010| ; IWilkins et al.l 1201 ll ) , when the universe was 
only about 500 Myr old. However, because most of these candidates are very faint (AB > 
27 mag), it is difficult to measure their spectra or measure their rest-frame optical SEDs to 
better constrain their stellar population and mass. It is imperative to search for relatively 
brighter (AB < 25 mag) galaxies at z > 6.5, and study them in great detail. There are three 
principal ways to accomplish this goal. First ap proach is through big cluster survey s, which 



are searching for lensed galaxies at z > 6.5 (e.g.. lHall et al.ll2012l ; iBradley et al.ll2012l ). These 
galaxies are intrinsically faint but because of the magnification their observed magnitudes 
are r elatively bright. Second approach is to take advantage of the HST parallel observations 



(e.g,. iTrenti et al.ll2011c lYan et al.ll201ll ) to explore large number of random fields at relatively 
shallow depth, and the final approach is to perform wide a r ea ground-based s urveys (e.g., 
Ouchi et all bood ; ICastellano et al.l bold : Irlickev et all bold ICapak et all boilh to identify 
and study such bright candidates. These methods are starting to identify bright galaxy 
candidates at z > 6.5 and are putting important constraints on the bright end of the rest- 
frame UV luminosity functions, though the number of these candidates remains too small 
to make any statistically significant conclusions about their physical properties. 

The wide area ground-based surveys have few distinct advantages compared to space- 
based approaches of identifying bright high redshift (z > 6.5) galaxies. First, they can cover 
much larger area than space-based observations, and secondly, the -fT s -band data, which is 
important for accessing the reliability of high redshift candidates, is only available through 
ground-based telescopes. These reasons make NIR ground-based surveys very appealing for 
searching bright high redshift galaxies. 
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Here, we present NIR data in the GOODS-N field obtained using the CFHT / WIRCam 
and the Subaru/MOIRCS instruments. We combine the archival data from these instruments 
to generate deeper images in J- and -fT s -bands. The Subaru/MOIRCS data cover ~60% of the 
central GOODS-N region in J, H, fQ-filters, while the CFHT/WIRCam covers a much larger 
area, extending well beyond the GOODS-N field in the J, i^-filters. The combination of the 
MOIRCS and WIRCam NIR data covers full area of the GOODS-N field (-169 arcmin 2 ), 
and is deeper in the central ~60% of the field. These combined images will extend far beyond 
the traditional field which the HST /WFC3 data will cover, and will considerably help to 
undertake multi-wavelength science goals that require NIR observations. As one application, 
here we use these NIR images along with the HST / ACS -B435, V606, ^775, ^850 data to identify 
bright Zs50-dropouts (i.e., LBG candidates at z > 6.5), and estimate their surface density in 
the GOODS-N field at the brighter limits (AB < 25 mag). 

We describe both the MOIRCS and WIRCam NIR observations along with existing 
HST / ACS data in § [21 the NIR data reduction is discussed in § [3], and the quality of our 
reduced NIR data is showed in § HJ We also discuss the -fT s -selected catalogs, their number 
counts and the completeness estimates in § HI In § El we develop and implement z 85 o-dropout 
selection criteria to identify LBG candidates at z > 6.5 and discuss the implications of our 
results. We summarize our results in § [6], emphasizing that we combine two NIR surveys 
of GOODS-N field to construct deep i^ s -band images for this field, needed for a variety of 
studies. 

Throughout this paper we refer to the HST /ACS F435W, F606W, F775W, F850LP fil- 
ters as -B435, V(306, ^775, ^850, and to the Spitzer /IRAC 3.6 jum, 4.5 //m, 5.8 //m, 8.0 jttm filters 
as [3.6], [4.5], [5.8], [8.0], respectively, for convenience. We assume a Wilkinson Microwave 
Anisotropy Probe (WMAP) cosmology with ^ =0.274, f) A =0- 7 26 an d H =70.5 km s" 1 Mpc" 1 



in accord with the 5 year WMAP estimates of lKomatsu et al.l (120091 ) . This cor responds to a 



look- back time of 12.93 Gyr at z ~ 7. Magnitudes are given in the AB system (jOke fc Gunn 



19831 ). 



2. Observations 

2.1. Subaru/MOIRCS Observations 

The MOIRCS is on the Subaru 8.2m telescope, providing wide-field imaging and spec- 
troscopic capability in the NIR bands. In the imaging mode, the MOIRCS provides 4x7 
arcmin 2 field of view with 0.117" pixel" 1 using two 2048x2048 HgCdTe (HAWAII2) ar- 
rays, each of which covers a 4x3.5 arcmin 2 field. There is almost no gap between FOVs 
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of each channel. Details of this instrument are given in llchikawa et al.l (120061 ). We used 
the Subaru archive SMOKA0 to retrieve the MOIRCS observations in the GOODS-N field. 
Typical exposure time for individual frames varied between ~50 to ~100 sees and standar d 



nine-point dither patterns (dithering by 10-15") were used (for details see, e.g.. iBundy et al. 



20091 ; iKajisawa et al.ll2011). These observat i ons were carried out over a period of thr ee years 



between 2005-2007 (e.g.. I Wang et al.l 120091 : iBundv et all 120091 : iKaiisawa et al.ll201lh . They 
cover an area of about 109 arcmin 2 , which overlaps with ~60% of the central region of 
the GOODS-N ACS field, as shown in Figure [TJ The typical seeing for these observations 
is between 0.6" and 0.8". Figure [1] shows the GOODS-N ACS field in gray shaded region 
with black solid line border, and the MOIRCS covered field with dashed black line. The full 
MOIRCS field is covered by 4 pointings and the field layout is shown by 2x4 chips (Figured]). 
These observations were done in three broad-band filters J, H and K s , whose transmission 
curves are shown as solid curves in the top panel of the Figure [2J 



2.2. CFHT / WIRCam Observations 



The WIRCam is on the CFHT 3.6m telescope and contains 4 2048x2048 pixel HAWAII2- 
RG detectors, covering a total area of 20x20 arcmin 2 with a sampling of 0.3" pixel -1 . 
WIRCam is a near-infr ared instrument sp anning the range 0.9-2.4 microns. Details of this 
instrument are given in IPueet et al.l 120041 ). We used the CFHT archive CADC@ to retrieve 
the WIRCam observations in the GOODS-N field. These observations were carried out over 



a period of three years between 2006-2008 (e.g., IWang et al.l 120101 ). They cover an area of 
about 1040 arcmin 2 , and overlap completely with the GOODS-N ACS field as shown by a 
large dot-dash box in Figured! The typical seeing for these observations is between 0.6" and 
0.8". These observations were done in two broad-band filters J and K s , whose transmission 
curves are shown as dash curves in the top panel of the Figure [2j Table [1] lists relevant 
information about the NIR data in the GOODS-N field, used in this paper. 

Filter transmission curves and total system throughputs used for these observations are 
shown in Figure [2j WIRCam system throughputs are calculated using their transmission 
coefficients^, and MOIRCS system throughputs ar e estimated by scalin g their transmission 
curves to the peak throughput given in Table 1 of llchikawa et al.l (120061 ). 



1 http://smoka. nao.ac.jp/ 

2 http:/ /cadcwww.dao.nrc.ca/cadc/ 

3 http:/ /www. cfht.hawaii.edu/Instruments/Imaging/WIRCam/WIRCamThroughput. html 
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2.3. HST/ ACS Observations 



The GO ODS-North field was ob served with HST /ACS broad-band filters -B435, V606, 
z 7 75, and ^850 (IGiavalisco et al.ll2004al ). The ACS imaging area covers about 169 arcmin 2 as 



shown in Figured! We used the GOODS team reduction of the ACS images (IGiavalisco et al. 
2004al ). These images have been drizzled from the original ACS pixel scale of 0.05" on to 



a grid of 0.03" pixels. We used v2.0 of the ACS GOODS images^ and corresponding AB 
magnitude zeropoints. The combination of deep HST /ACS S435, Veo6, Z775, z$5o and ground- 
based J, K s observations provide the wavelength coverage needed for selecting LBGs at 

> 



z £ 6.5 (Figure [2]). 



3. Near-IR Data Reduction 



We have reduced publicly available WIRCam and MOIRCS NIR imaging data in the 
GOODS-N field by using a reduction pipeline based on the Interactive Data Language (IDL) 
- Simple Imaging and Mosaicking PipeLinE (SIMPLE) developed by Wei-Hao W ang. Details 



of this pipeline and the steps required to reduce the WIRCam data are given in IWang et al. 



(12010!) . We have independently reduced the WIRCam data using this pipeline. Reducing 
the MOIRCS data is less straightforward and we took the following steps, using the SIMPLE 
software. Many general steps are similar for both datasets. 



• Generate File lists: We selected all archival observations taken on photometric nights 
only. SIMPLE requires input lists of the raw FITS files which are grouped in such a 
way that each group has same exposure time, is from the same chip, is taken on the 
same day, and has same coordinates. If a large group of exposures is taken at the same 
pointing then we break them up into smaller chunks. 

• Reduction: For the first pass reduction, these file lists are given as input to the main 
IDL procedure reduce_moircs .pro included in the SIMPLE software. This script does 
basic reduction steps including the flat-fielding, cosmic ray removal, sky-subtraction, 
aligning images, distortion corrections, and general plate solution. It also produces 
exposure time maps. 

• De-fringing: The output of the first pass reduction is carefully checked for image defects 
and/or fringes. The MOIRCS data show strong circular fringes in both chips but more 



4 http://archive. stsci.edu/pub/hlsp/goods/v2/ 
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prominently in chip-2. We use the program def ringe_moircs .pro provided in the 
SIMPLE package to remove the fringes. The program first masks detected objects in 
a flattened and sky-subtracted image and transforms the image to a polar coordinate 
system, where the circular fringes become nearly straight lines. Then it fits a 5th 
degree polynomial to each column along the fringes to produce a model fringe image. 
The fringe model is then transformed back to the original Cartesian coordinate and 
subtracted from the image. This procedure considerably improves the quality of the 
reduced images. A very similar fringe r emoval procedure is adopted by the MOIRCS 



Deep Survey team (IKajisawa et al.ll201l[ ). A few exposures with very bad fringes, which 
could not be removed through above mentioned process, and/or image defects are not 
included in the second pass reduction. The reduced images and the exposure time 
maps are used to form large mosaics. 

Mosaics: The final combined mosaic is generated by the procedure mosaic_wide .pro. 
The final MOIRCS images in three filters are registered to each other and to the 
GOODS ACS images by making the WCS of the MOIRCS images consistent with 
other multi-wavelength images of the GOODS survey. The final pixel scale for the 
MOIRCS images is 0.15" pixel -1 . For proper comparison with the WIRCam 0.30" 
pixel -1 images, we also generated 0.30" pixel -1 MOIRCS images using the IRAF task 
'blkavg'. 

Zeropoints: The procedure reduce_std_moircs .pro was used along with the standard 
star images from the archive to reduce standard star images and calculate zeropoints 
in all the three MOIRCS filters (J, H, K s ). The MOIRCS images are in ADU/sec and 
their derived AB zero-points are 26.175 mag for K s , 26.569 mag for H and 25.965 mag 
for J filters. The WIRCam images are in /iJy units with a zeropoint of 23.9 mag. 

Combined Images: We have generated combined (WIRCam and MOIRCS) J, K s im- 
ages, whi ch have similar see ing/ point spread functions (PSFs), using the Terapix SWarp 



package (IBertin et al.ll2002l ). which resamples and co-adds FITS images using astromet- 
ric projection defined in the WCS. The goal of these combined images is to generate the 
deep J and -fT s -images in the central ~60% of the GOODS-N field. The combination 
was performed as follows: First the images were scaled to a common flux scale. Abso- 
lute root mean squared (RMS) noise maps were generated for each mosaic by scaling 
the mean variance of the images to the mean inverse exposure time maps generated by 
the SIMPLE reduction package. The data were then combined using a weighted mean 
combination with the SWarp package using the RMS images as the weight image. A 
bi-linear interpolation kernel was used to re-sample the images onto a common pixel 
grid. The final combined J, and K s images are in /iJy units and have a zeropoint of 
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23.9 mag. 



4. Near-IR Data Quality and Catalogs 

We have reduced archival J, _K" s -band WIRCam and MOIRCS images, and shallow H- 
band MOIRCS image of the GOODS-N field. We then generated combined images (WIRCam 
and MOIRCS) in the J and i^ s -filters, which covers much larger area than the HST data. 
Table [1] shows area coverage information of the NIR data in the GOODS-N field. In this sec- 
tion, we discuss the quality of these combined images, source detections, and the generation 
of catalogs covering the entire area of this field. 



4.1. Astrometry 

The astrometric solutions of the final mosaiced images were tested by comparing them 
with the available HST /ACS zsso-band catalogs. Figure [3] shows the relative astrometric 
offsets in RA and DEC between the _ff s -selected and the GOODS-N ACS z 850 -selected cata- 
logs. For the purpose of this comparison, we selected well detected/non-saturated compact 
sources with S/N > 20 and full width half maximum (FWHM) < 1.2" from the GOODS-N 
ACS ^so-selected catalog^ and matched them with the if s -selected catalog from the com- 
bined image. Figure [3] shows the distribution of the offsets with their mean, median and 
sigma values. The uncertainty (a) in the offsets between the combined _K" s -selected and the 
ACS 2; 850 -band catalogs is 0.19" for RA (uncertainty in RA-cos(DEC) is 0.09) and 0.09" for 
DEC. The uncertainty in the RA is slightly higher than the DEC, but these are still less 
than one pixel in the i^ s -band (0.30" pixel -1 ). These uncertainties include both the internal 
errors in the ACS ^gso-band catalog and the errors introduced during the image registration 
of SIMPLE. There are no systematic offsets between the K s positions and the ACS 
catalog. 



4.2. Photometry 

To check the consistency of the photometry in our reduced images, we compared the K B - 



band magnitudes obtained using the SExtractor (IBertin fc Arnoutdll996l ) AUTO apertures 



in the WIRCam image with the MOIRCS image, and also compared i^ s -band magnitudes 



5 http: / / archive.stsci.edu/pub/hlsp / goods / catalog_r2 / 
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in the combined image with the K s -band magnitudes from the Two Micron All Sky Survey 
(2MASS). The 2MASS magnitudes have b een converted to the AB magnitude system using 
the conversion given by lCiliegi et al.l (120051 ) . The top panel of Figure H] shows the comparison 
between the SExtractor K s magnitudes, for all matched objects, from the WIRCam and the 
MOIRCS images. There is no significant systematic offset between these two magnitudes, 
and the mean of their difference is ~0 mag, with the uncertainty in the difference being 
-0.30 mag. The mean value (uncertainty) for [K s (WlRCam)-K s (MOIRCS)] is 0.03 (0.08) 
for brighter objects (K s <20 mag). The bottom panel of Figure H] shows the comparison 
between the SExtractor K s magnit udes, measured i n the combined images, and the K s 
magnitudes from the 2MASS catalog (jCutri et al.ll2003l ). We compare these magnitudes over 
the range ~14 to ~16 mag (as shown by dot-dash vertical lines on the bottom panel of the 
Figure H]) because of the non-linearity issu es in the WIRCa m and the selection effects in 
2MASS outside this magnitude range (see IWang et al.l |2010| . for details). The comparison 
here is based on a small number of bright sources and shows no major systematic offset, 
with the uncertainty in the difference being very small and mean magnitude difference of 
0.02 mag. 



4.3. Catalogs and Number Counts 

We used the SExtractor in dual-image mode with the K s image as the detection im- 
age for /^-selected catalog, and J image for J-selected catalog. We generated the RMS 
maps from the weight /exposure time maps and used them during the SExtractor runs to 
get a better estimate of the photometric errors, and hence, of the S/N. Here, S/N is de- 
fined as 1.0857 divided by the SExtractor error in the total magnitude, which is similar to 
FLUX_AUTO/FLUXERR_AUTO. These RMS maps also help to exclude false and spurious 
detections at the edges with low signal-to-noise values. We set the DEBLEND_MINCONT 
parameter to 0.0001, DETECT_MINAREA to 4 pixels (2 pixels for J-selected catalog) and 
detection threshold to 2. Oct (la for J-selected) , and used convolution with a 3 pixel Gaussian 
filter in the SExtractor configuration file. Table [2] lists the SExtractor parameters used 
to generate K s and J-band selected catalogs in the GOODS-N field. We used SExtractor 
MAG_AUTO as a measure of the total magnitude. The left 4-panels in Figure [5] show the 
source counts (i.e., number per arcmin 2 per 0.5 mag bin) in K s and J bands. The leftmost 
two panels show the MOIRCS and WIRCam /A^-band counts. They both reach comparable 
depths marked by dashed vertical lines. The right panels show the K s - and J-band number 
counts obtained from the combined image. Because of the non-uniform exposure time in the 
WIRCam and MOIRCS images, the combined i^ s -band image is ~0.3-0.5 mag deeper than 
either WIRCam or MOIRCS image. 
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We quantify the depth and completeness in our K s - and J-band images by inserting 
numerous fake sources of varying magnitudes into the reduced images and recovering them 
using the same SExtractor parameters (Table [2]) used for the real sources. Fake sources 
were given Gaussian profiles with the FWHM values similar to the point sources in our 
image. We define the image depth by the magnitude corresponding to a recovery rate of 
50%. This depth implies a slightly lower detection rate for real galaxies, which are more 
extended. The rightmost two panels in Figure [5] show the 3a completeness curves (dashed 
curves) for both images. The top panel shows completeness for .K" s -band image, while the 
bottom panel is for the J-band. The dashed vertical line in each panel shows the magnitude 
at which the number counts fall to 50% of their peak value, which we consider as a measure 
of the completeness limit for these images. 

We also quantify reliability of our catalogs by doing the negative image test. We multiply 
the science (positive) image by -1 to get a negative image. We then apply same detection 
procedure using SExtractor parameters listed in Table [2j The number of detected 'objects' 
in these negative images quantify reliability of these catalogs. The dot-dash curve in the 
rightmost two panels of Figure [5] shows the reliability curve for both J- and Ag-selected 
catalogs. The dot-dash vertical line shows the magnitude at which the reliability falls to 
50%. 

The simulations and the observed number counts show that both (J and K s ) combined 
images are complete to AB magnitude limit of ~25 mag (3cr), and 50% reliable to ~24.5 mag. 
We also have if-band data in the GOODS-N region from the Subaru telescope observations. 
The if -band MOIRCS image is shallower compared to J or A" s -bands by ~1 mag, but can 
be very useful for brighter objects in this field. 



5. Lyman Break Galaxies at z > 6.5 (z 850 -dropouts) 

As an application of these NIR data in the GOODS-N field, we combine these with 
the existing HST /ACS data to look for z 850 -dropout galaxies, which are LBG candidates 
at z > 6.5. To identify ^so-dropout s, we search for Lyman-break signature that occurs at 



the rest-frame 1216 A (IMadaul 119951 ) using the dropout technique. This technique is based 
on photometric color selection and requires imaging in at least two filters, one to the blue 
side of the break and the other to the red side (usually we have two or more filters on both 
sides). The presence of Lyman break makes high-redshift galaxies much fainter in the blue 
band than in the red one, or in other words, it makes them appear to "drop-out" from 
the blue band. For this reason, this method is also known as the dropout selection, and 
the candidates found in this way are generally referred to as "dropouts". At z > 6.5, this 
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break moves through observed ^850-band (~9100 A). Figure [2] shows the locations of the 
rest-frame 1216 A Lyman break at z > 6.5. It is clear that ground-based J and K s band 
filters, combined with the ACS bands are crucial in identifying LBG candidates at z > 6.5. 

The drizzled ACS images in -B435, Vm§} 2775, ^850 filters have a spatial resolution of 0.03" 
pixel -1 , compared to combined J and A" s -band images which have 0.30" pixel -1 resolution. 
Therefore, to perform matched aperture photometry, we generated re-binned (10x10), PSF 
matched, ACS images to properly compare with combined J and -fr s -band images. Details 
of catalog generation, completeness, and reliability are discussed in § 14.31 In summary, we 
performed matched-aperture photometry on -B435, V606, ^775, ^850, J, K s images by running 
the SExtractor in the dual-image mode with J-band as the detection image. We used 
corresponding RMS maps to reduce the number of spurious sources detected in the low 
signal-to-noise ratio edge of the images. For object identification, we adopted a limit of 
at least 2 contiguous pixels above a threshold of la. We constructed a J-selected catalog 
using SExtractor parameters shown in Table[2j We used SExtractor MAG_AUTO, with the 
default Kron factor=2.5 and minimum radius=3.5, as a measure of the total magnitude. This 
catalog was used to select LBG candidates at z > 6.5 based on the dropout color selection 
criteria. The colors of the objects were estimated from the matched-aperture MAG_AUTO 
photometry on the HST /ACS and combined NIR images, as discussed above. 



5.1. Zg5o-dropout Selection 



5.1.1. Color Selection 



The ^850- drop out selection uses color criteria obtained from the stellar population mod- 
els of iBruzual h Charlotl (120031 . BC03). Figure [6] (top panel) shows the BC03 star-forming 
galaxy models with E(B-V)=0, 0.15, 0.30 mag (solid blue lines with open circle s show- 
i ng di fferent redshifts), expected colors of late-type stars (black filled circles) from Pickles! 
( 19981). and tracks of lower- redshift ellip ticals (red lines) from three different models (BC03, 
Kinney et al.l Il996l ; IColeman et al.l 119801 ). Our selection region is shown as a gray shaded 
area. We have corrected the colors for IGM attenuation using the prescription of iMadau 



( 119951 ) . The two-color dropout selection criteria ado pted here are similar t o those used to 



ident i fy and study LBG candidates at z g 1-8 (e.g ., 



Giavalisco et al. 2004b; Bouwens et al. 



20081 : lHathi et al.ll2008l : lYan et al.ll20ld : lHathi et aUBOlOh . The J-selected catalog is used to 
select zsso-dropouts from the (zgso _ J) vs. (J - K s ) color-color diagram (Figure [6]) using 
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the following criteria: 



Os50 - J) > 1.5 mag 

and (J — K s ) < 1.2 mag 

and (2 850 - J) > 0.99 + [0.85 x (J 

and (J — K s ) > —0.5 mag 

and S/N(J) > 2.0 



mag 



We also require non-detection of the candidates (S/N<1) in ACS optical -B435, V606, 2775- 
bands. The S/N cut in the J-band (> 2) corresponds to an AB magnitude limit of ~25.3 mag. 
We have also applied the additional criterion (J - K s ) > -0.5 mag to eliminate the possibility 
of selecting spurious candidates (e.g., ex treme equivalent width objects at lower redshifts: 



van der Wei et al.ll201ll ; lAtek et al.ll201ll ). since this color-color space is based only on z 850 , 
J, and K s filters, and it is required that LBG candidates be detected and are relatively 
bright in the i^ s -band. After initial color and S/N cuts, we visually inspected each selected 
object to remove any spurious source (due to image defects or their proximity to a bright 
foreground object, edge effects, or faint stellar diffraction spike). This leaves 14 potential 
candidates as shown in the top panel of the Figure [6j 



5.1.2. Redshift Selection Functions 

To generate redshift selection functions, we ran simulations to calculate P(m, z), which 
is the probability that a galaxy with apparent magnitude m at redshift z will be detected 
in the image and will meet our color selection criteria (§ 15.1.11) . In these simulations, large 
numbers of artificial objects with a range of redshifts and magnitudes were added to real 
images, and then recovered using exactly the same method and selection criteria that were 
employed for the real observations. We simulated these objects in four HST / ACS, J, and 
K s bands because we use these filters to perform dropout color selection as shown in § 15.1.11 
The input spectrum for each simulated object is from BC03 models with constant star- 
formation his tory, three different E(B-V) values, two different metallicities, and applying 



Madaul (119951 ) prescription for IGM attenuation below rest 1216 A. The bottom panel of the 



Figure El shows the redshift selection function for different J magnitudes. 



5.1.3. Spitzer/7iM C Colors 



The GOODS-N field has deep Spitzer/IRAC imaging in [3.6], [4.5], [5.8], and [8.0] chan- 
nels. These observations are vital in assessing the reliability of our ^o-dropout candidates. 
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The (J - [3. 6]) and {J - [4. 5] ) colors are used to differentiate between extremely red ob- 
jects (EROs; lYan et al.l l2004h at z ~ 2 and high redshift candid ates. The EROs have (J 
[3.6]) and (J - [4.5]) colors in excess of ~2.5 (jYan et al.l 120041 ) . putti ng a constraint o n 



the colors of high redshift galaxy candidates to have bluer colors (e.g., lEyles et al.l 120071 ). 
To examine this, we used Spitzer/IRAC [3.6] and [4.5] images to estimate magnitudes of 14 
candidates in these bands. We simply calculated IRAC total magnitud e s from 3"-diameter 
aperture magnitudes and used aperture corrections given by lYan et al.l (120051 ) . To confirm 
our photometry across different NI R passbands, we com pared our SExtractor magnitudes 
with the Template FITting (TFIT; lLaidler et al.l 120071 ) magnitudes. The TFIT technique 
performs consistent multi-waveband photometry on images with widely different resolutions. 
Our J, K s , [3.6], and [4.5] SExtractor magnitudes are very similar (within uncertainties) 
to the TFIT magnitudes obtained using different set of proprietary CFHT NIR images in 
this region. To confirm our color selection and to differentiate between LBG candidates at 
z > 6.5 and EROs at z ~ 2, we estimated photometric redshifts for 14 candidates. 



We used three different codes t o estimate photom etric redshifts, name l y HyperZ (|Bolzonella et al 



20001 ). EAZY (iBrammer et all 12008 ). and Le PHARBlj (lArnouts et al.lll999l : lllbert et al.l 120061 ). 
Out of 14 candidates, the photometric redshift estimate for 12 objects was not consistent 
between three codes with few objects showing possible detection in optical bands during 
visual inspection. Therefore, we find two possible LBG candidates at z > 6.5 after applying 
our color criteria, visually confirming non-detections in optical bands, and based on our pho- 
tometric redshift estimates. These two candidates are shown in Figure Table E] lists the 
coordinates and TFIT photometry of these candidates. Their spectral energy distributions 
and best-fit models based on Le PHARE code are shown in Figure El Though photometric 
redshifts show high probability that these are LBG candidates at z > 6.5, the exact nature 
of these LBG candidates is unclear, and we will explore their properties and few possible 



scenarios. 



5.2. Are these luminous objects at z > 6.5? 



Based on the photometric data, it is unclear whether the objects found here are genuinely 
at z > 6.5 or not. Both objects meet th e color selections suggested by previous studies (e.g., 



Bouwens et al.ll2008l ; ICapak et al.ll201ll ). but they are very bright ( J < 24 mag) to be at this 
redshift, specially in a comparatively smaller area (~169 arcmin 2 ) of the GOODS-N field. 
In this section, we discuss properties of these candidates. 



http://www.c£ht.hawaii.edu/^arnouts/LEPHARE/lephare.html 
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5.2.1. Morphologies 

The GOODS-N catalog of z 7 7 5 -band dropouts^ selected using criteria outlined in Beckwith et al 



( 120061 ) includes these two candidates, but this catalog is only based on HST /ACS -B435, Veoe, 
^775) ^850 magnitudes. Addition of J and K s observations help to distinguish between an 
object at z ~ 6 and z > 6.5, depending on the location of Lyman break in the Zgso-band. 
Both these candidates are very faint but compact in z 8 5o-band (stellarity > 0.5) so there is 
a strong possibility that these objects are faint point sources (late- type stars). The J-band 
SExtractor stellarity for candidate A is 0.54, while it is 0.13 for candidate B, but the J- 
band PSF is less sharper than ACS so J-band stellarity could be a less convincing measure 
of compactness for these objects. Therefore, based on z 850 -band compactness measurements, 
it is possible that these zgso-dropout candidates could be stellar-like. 



5.2.2. Passive Galaxies 

One of the major interlopers in selecting high redshift dropout galaxies is the lower 
redshift passively evolving galaxies. The 4000 A Balmer break in old galaxies at z ~ 2 
could potentially be mistaken for a Lyman break at z > 6.5. These passive galaxies can be 
modelled with suitable set of spectral synthesis models. We use the BC03 models to predict 
the colors of such objects. Figure [6] (top panel) shows predicted colors of these old galaxies as 
a function of redshift through various evolutionary tracks plotted in red. The color selection 
region of z 850 -dropouts (gray shaded region in the top panel of Figure based on (z 850 - 
J) and (J - K s ) avoids these tracks, and hence, minimizes contamination of our sample by 
these passive galaxies. Also, the EROs at z ~ 2 have very red (J - [3.6]) and (J - [4.5]) 



colors in excess of ~2.5 ( lYan et al.ll2004j ) while our z 850 -dropout candidates have much bluer 



NIR colors, which also suggests that these candidates are not likely EROs at z ~ 2. 



5.2.3. Late-type Stars 

The NIR colors of Galactic cool stars are also very similar to the high redshift dropout 
candidates and can mimic the Lyman break in the color-color space as they are also optically 
faint and bright in the NIR. The main distinguishing properties between a high redshift 
galaxy candidate and a cool dwarf are their morphologies and NIR colors. We discussed 
morphologies in § 15.2.11 Here, we investigate the J, K s and Spitzer /IRAC colors which will 



7 http:/ /archive. stsci.edu/pub/hlsp/dropouts/ 
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shed light on the nature of these zgso- drop out candidates. The (zs5o - J) and (J - K s ) colors 



i n Fig ure [6] show expected colors of late-type stars (black filled circles) based on iPickles 



( 119981 ) models. The colors of candidate B fall in the region occupied by the cool dwarfs, 
while candidate A has redder (J - K s ) color. Therefore, based on these colors, it is possible 
that candidate B could be a late-type star. 



Stanway et al.l (120081 ) have proposed J, K s and Spitzer /IRAC color cuts, which were 
designed to remove passive galaxies at z < 2 and late-type stars. The 2 850 -dropout candidate 
A successfully passes these four color criteria, while candidate B does not pass the ([3.6] 
[4.5]) color cut. These criteria also indicate that the candidate B could be an interloper 
(most likely a late-type star) based on these NIR colors. 



5.2.4- Photometric Redshifts 



Kajisawa et al.l ( 120111 ) have analysed MOIRCS data, separating it into two regions — the 
deep and the wide. Our zgso-dropout candidat es are not locate d in th e deep MOIRCS region, 
but are covered by the shallower wide region. Kajisawa et al.l ( 1201 ll ) estimated photometric 
redshifts of these candidates using their own MOIRCS wide photometric catalog. They 
measured photometric redshift using four different techniques. Three out of four techniques 
estimate redshift of these candidates to be at z > 6.5, while the fourth approach (using EAZY 
code) estimates lower redshift (z — 2) for these candidates. 

We performed SED fitting using available multi-wavelength data (optical to NIR) in the 
GOODS-N field to predict photometric redshifts of our z 850 -dropout candidates. We used 
three photometric redshift codes (Le PHARE, HyperZ and EAZY) on TFIT photometry to 
estimate t heir redshifts . The TFIT photometry generates accurate photometric redshifts, as 
shown by iDahlen et al.l ( 120101 ) for the GOODS-S field. These codes predict high probability 
for these two candidates to be at z > 6.5. Figure [8] shows SEP fitting using Le PHARE code 
for two types of templates — blue for a late-type star (from lChabrier et al.ll2000l library) and 
red for a star forming galaxy (from BC03 library) at z > 6.5. The SED fitting for candidate 
B shows equal probability (very similar \ 2 ) f° r both these templates, while candidate A has 
a higher probability for a LBG at z > 6.5. This implies that candidate B could be a late- type 
star, but it is also equally likely to be a high redshift galaxy. We also note that due to lack of 
y-band data, the exact location of the Lyman break between the £350- and J-band could be 
somewhat uncertain, and that would increase the uncertainty in the photometric redshift. 
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5.2.5. Keck Spectroscopy 



We used the DEep Imaging Multi-Object Spectrograph (DEIMOS; iFaber et all 12003 ) 
on the 10m Keck telescope to perform spectroscopic observations of these two candidates 
during our March-April 2011 observing run. The total exposure time for these observations 
was ~4 hours. The seeing was in the range of 0.5" - 0.7". We used the GC455 filter and 
the 600 line mm" 1 grating, blazed at 7400 A. The spectral coverage was between ~5500 A 
and ~10000 A. The spatial pixel scale was 0.1185" pix -1 , and the spectral dispersion was 
0.65 A pix -1 . The slit widths were 1". The preliminary r eduction were performed using the 



spec 2d IDL pipeline^] developed by the DEEP2 Team ( Cooper et al. 2012 ; Newman et al. 



20121 ). Wavelength calibration was performed by fitting to the arc lamp emission lines. The 
standard stars and science targets were observed using the same slit width so we did not 
apply any corrections due to slit-loss effects. 

We present the Keck spectrum of one z 8 50-dropout candidate which has high enough 
S/N ratio to make a redshift measurement possible. For the second candidate, because of 
the unfortunate location on the mask and sky lines, the spectrum is not reliable. In Figure [9j 
we show one and two dimensional DEIMOS spectra of the ^sso- drop out candidate B between 
8400 and 10000 A. Figure |9] shows a weak emission line at ~9090 A indicated by a dot-dash 
line in the one dimensional spectrum, and by a red arrow in the two dimensional stamp. 
The sky spectrum is shown at the bottom of the one dimensional spectrum and there is no 
major sky line at this wavelength. This spectral feature also shows asymmetric nature of 
a Lya line, and we don't see any clear detection of other lines, so we tentatively identify 
this line as a Lya at z ~ 6.5. The estimated line flux is ~1.4xl0 -17 erg s _1 cm -2 , and 
the rest-frame equivalent width (EW) is >21 A. Hence, candidate B is a possible LBG at 
z ~ 6.5, consistent with its photometric redshift. The red part of the spectrum of candidate 
A is mostly affected by improper masking and sky lines, so it is unreliable and inconclusive. 
The bluer part of the spectrum (below 8000 A) does not show any clear detection of emission 
line. This could imply that candidate A is not a low redshift galaxy, but does not rule out 
other possibilities (e.g., late- type star or LBG). This spectrum cannot conclusively predict 
the nature of candidate A — whose photometric redshift is ~7.4 — because the possible 
Lya emission line would be around 10200 A, which we cannot identify. Therefore, the 
Keck spectrum of candidate A cannot differentiate between a late-type star or a LBG at 
z > 6.5. We will try to re-observe these candidates with Keck/DEIMOS or MOSFIRE NIR 
spectroscopy. 



The three spectroscopically confirmed objects of lOno et al.l (120121 ) at z > 6.5 were also 



3 http: / /astro. berkeley.edu/cooper / deep / spec2d/ 
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observed during the same Keck observing run. These objects were about 1 mag fainter in 
the continuum than our candidates, and all three showed detectable Lya in emission. These 
successful Lya detections along with the likely Lya detection in Candidate B allo w us to put 



l ower limits on Lya line emission from our z 85 o-dropout candidate A. The three lOno et al 



( 120121 ) objects have Lya flux limits (la) in range of 2.5-2.7x10 17 erg s 1 cm 2 , and rest- 



frame EW between 33 and 43 A. These flux and EW measurements are consistent with other 



confirmed Lya emitters at z > 6.5 (e.g., IVanzella et al.l 1201 ll ; iPentericci et al.1 1201 ll ). This 
implies that if our candidate A is a Lya emitter, then we could detect Lya line with flux 
> 1.4xl0~ 17 erg s" 1 cm" 2 , and EW>20 A. 

The fraction of LBGs showing Lya in emission is a strong function of the absolute UV- 



magnitude (Muv) 
observations at z 



i.e., fraction increases for 

i 6-8 



(e.g., 



Stark et al. 2011 



particle hydrodynamics simulations (e.g., iForero- Romero et al. 



ainter galaxies. This has been confirmed by 



Schenker etal 



2012h. and also by smoothed 



20121 ) . These studies imply 



that more luminous LBGs are less likely to show Lya in emission, and it is possible that 
candidate A will not show Lya line bec ause of its luminos ity. We should also emphasize 
that, at similar continuum magnitudes, ICapak et al.l ( 120111 ) could not clearly detect Lya 
emission for their LBG candidates at z > 6.5. Therefore, to clearly identify the nature of 
our LBG candidates at z > 6.5, we need deep NIR spectroscopy from Keck/MOSFIRE or 
HST/WFC3. 

In this section, we investigated morphology, NIR/IRAC colors, and Keck spectroscopy 
of our two 2 850 -dropout candidates. We find that though photometric redshift and NIR 
colors of candidate A support that it could be a LBG at z ~ 7, because of inconclusive 
Keck spectrum the true nature of this candidate is uncertain. The candidate B could be a 
real Zg5o-dropout galaxy at z ~ 6.5 based on a weak Lya detection in the Keck spectrum, 
but we cannot rule out other possibilities based on our data. Therefore, it is plausible that 
one or none of these candidates are at z > 6.5. Hence, we conclude that we find 1±1 LBG 
candidates at z > 6.5, brighter than AB~24.5 mag, in the GOODS-N field. 



5.3. Discussion 



5.3.1. Other Surveys and Cosmic Variance 

The rarity of bright LBG candidates at z > 6.5 in a comparativel y smaller area 
(~169 arcmin 2 ) questions as why other surveys have not found such objects. IConselice et al. 



( 1201 it ) used 30 NICMOS pointings of massive galaxies at 2 ~ 1.7-2.9 to cover smaller 
area of the GOODS-N field with NIR data. These NICMOS observations were used by 
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Bouwens et all (l2010bl ) to search for LBGs at z > 6.5, but th ey do not cove r the area in 
which we find our two candidates. The PEARS grism survey (jMalhotral 120071 ) in this field 
also does not cover this region, which could have helped to confirm the nature of our z^q- 
dropout candidates. 



Ouchi et al.l (120091 ) survey was a F-band survey in the GOODS-N field and does not 
use the NIR bands. This survey requires a detection brighter than 26 mag in the F-band 
(~1 micron). Their redshift selection function shows a survey optimized to select galaxy 
candidates at 6.5 < z < 7.1, therefore if our candidates are outside this redshift range then 
Ouchi et al. J2009f ) could not have found them. Moreover, they used Suprime-Cam on the 
Subaru Telescope to do a wide-area survey at 1 micron, which is not possible in the NIR 
simply because of the small field of view of NIR detectors. As a result, they also surveyed the 
extended GOODS-N region with no HST data. In the present study, we have generated deep 
NIR images in the GOODS-N area for which we have very deep HST /ACS and Spitzer /IRAC 
data, all with consistent photometry using TFIT technique. Combining these data with the 
NIR data, we can then search for the ^850-band dropout galaxies, which puts limits on the 
source density of these sources. Also, the availability of combined -fT s -band data allows 
searches for high-z galaxies when HST /WFC3 (YJH ) data become available through the 
CANDELS jGrogin et al.lboill ; iKoekemoer et al-lboilh . 



Similar wide area NIR surveys have been done in the GOODS-S and other fields to look 
for brighter LBG candidates at z > 6.5. iHickey et al.l ( 120101 ) used ESO/VLT observations in 
the GOODS-S field covering ~119 arcmin 2 to search for high redshift LBGs. Th ey found four 



possib le candidates at z > 6.5 with an AB magnitude brighter than 25.5 mag. IHickey et al. 



( 120101 ) survey covers slightly smaller area than our survey and does not find any 2; 850 -dropouts 
with J magnitude br ighter than 24 mag, but they do find objects with similar K s magnitudes 
as our candidates. ICastellano et al.l ( 120101 ) also conducted a VLT survey covering three 
different fields (total area 161 arcmin 2 ) to search for LBGs at z > 6.5. They found a total of 
15 candidates but all are fainter t han 25.5 mag. The HST yWFC3 parallel imaging surveys 
( lYan et al.ll2011t iTrenti et al.ll201lf ) have covered ~125 arcmin 2 area — which is smaller than 
our current survey area — to search for Y 098 -dropouts, which is optimized to identify LBGs 
at z > 7.5. These surveys have identified few candidates with AB magnitude fainter than 
or equal to 25 mag. Though most of the surveys have smaller area coverage than our full 
area GOODS-N field (~169 arcmin 2 ), they don't find any LBG candidates at z > 6.5 with 
magnitudes brighter than Jab ~24.5 mag. 



Stanway et al.l ( 120081 ) performed extensive search for LBGs at z > 6.5 in 10 widely 



separated fields as a part of the ESO Remote Galaxy Survey (ERGS). This survey covered 
~225 arcmin 2 and found one possible LBG candidate at z > 6.5 with Jab ~23.6 mag, which 
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is very similar to J magnitude of our zgso-dropout candidates. If IStanway et al.l (120081 ) bright 
candidate is confirmed, then the surface density of such bright LBG candidates at z > 6.5 
is ~0.4xl0~ 2 per arcmin 2 in the ERGS. 

Based on other wide area surveys, it seems very unlikely that both of our bright candi- 
dates are truly at z > 6.5. The GOODS-N field is a comparatively sm all area to find two 
such bright LBG candidates at such a high redshift. ICapak et al.l ( 120111 ) found three possible 
z 850 -dropout candidates with similar brightnesses, but in a area as large as 2 deg 2 . That 
being said, if one of our candidates is confirmed to be at z > 6. 5 then the estimated surface 
density is ~0.5xl0 -2 per arcmin 2 , which is very similar to what IStanway et al.l (120081 ) found 
in the ERGS. Therefore, finding such a bright LBG candidate at z > 6.5 in a survey area of 
~169 arcmin 2 is a statistical possibility based on cosmic variance. 



5.3.2. Implications 



We find two possible LBG candidates at z > 6.5 from our analysis of NIR images in the 
GOODS-N field (~169 arcmin 2 ). The first candidate could be a real LBG at z ~ 6.5 based 
on a possible Lya line in the Keck spectrum, and the photometric redshift of the second 
candidate support its high redshift nature. If confirmed, these are among the brightest LBG 
candidates at z > 6.5. At a redshift of z ~ 7, we estimate absolute magnitudes Muv —-23.7 
and -23.1 mag for candidates A and B, respectively, which is ~1 mag brighter than spec- 
troscopically confirmed Lya emitting z 85 o-dropouts by lOno et al.l (120121 ) . The two possible 
LBG candidates at z > 6.5 in the GOODS-N field implies a surface density of ^l.OxlO" 2 



per arcmin 2 at J < 24.5 mag. This i s consistent with other surveys (e.g., IStanway et al. 



20081 ; lOuchi et al.l 12009b lYan et al.l 1201 11 ) focusing on galaxies brighter than 25.5 mag. 



Using iMadau et al.l (119981 ) formalism, the rest-frame UV luminosities of our candidates 
correspond to star formation rates (SFRs) of ~100-200 M Q yr _1 . As a rough estimate, if they 
keep forming stars at the same rates to z~5, they could accumulate stellar masses of (0.5- 
1.0) xlO 11 M in next -400 million years. Therefore, our luminous ^Bo-drop out candidates 



could be progenitors of massive LBGs (~10 n M ) observed at z~5 (e.g., IWiklind et al. 



20081 ). 



The bright-end of the luminosity function declines exponentially so we expect fewer 
bright galaxies at these redshifts. The surface density, which we predict based on our lumi- 
nous LBG candidates at z > 6.5, is significantly higher than what is expected from smaller 
area fainter HST surveys, but such a difference at the bright end of the lumino sity function is 



not surprising and is within the theoretical limits based on halo mass function ( IBouwens et al 
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20081 ; ICapak et al.ll201lf ). The star formation processes are not well understood at these red- 
shifts, so it is very likely that the mechanism involving truncation of star formation and 
feedback are different at higher redshifts and this can lead to a flatter bright-end for the 
luminosity function. 

To estimate the number density of these possible bright ^sso-dropouts, the GOODS-N 
volume was estimated using P(m, z) simulations (details in § !5.1.3p . This volume estimate is 
yery s imilar to the estimate obtained using cosmic variance calculator of iTrenti fc Stiavelli 
(120081 ) . The number density for one object at z > 6.5 obtained using t his volume is plotted 



in Figure ITUl along with other number densities from the lite rature ( IBouwens et al.l 12008 



Oesch et al.l I2009L BOld : ICapak et al.l 1201 ll ; Bowler et al.l 120121 ) . The number densities at 



M 



uv 



^-23. 5 mag, if confirmed, point to a 'no-knee' type luminosity functi on at z > 6.5, 



which has been consider ed as a possibility for z > 6 luminosity functions (e.g.. IBouwens et al 
20081 ; ICapak et al.ll201ll ). The bright-end of the luminosity function cannot be constrained by 
faint, small area surveys. To find brightest galaxies and get better constraints fo r the bright 



end of the LF at these redshifts we wil l need wide area, shallow surveys (e.g., ITrenti et al. 



2011; 


Yan et al. 


2011; 


Groein et al. 


2011) 



difficult to get spectroscopic confirmations — even with 10m Keck telescope ( as discussed 
in § I5.2.5P — if they are not strong Ly-a emitters (also see ICapak et al.l 1201 ll ). The high 



resol ution HST imaging in NIR with the CANDELS and deep WFC3 /IR grism observations 



(e.g., 



Brammer et al. 



20121 ) may help to confirm these and other such candidates. 



6. Summary 

We have combined archival NIR images obtained from the Subaru and the CFHT tele- 
scopes to generate deep J and K s band images in the GOODS-N field. These NIR data 
are vital for variety of multi-wavelength science goals. We use these images to search for 
bright z 8 5crdropouts, (i.e., LBG candidates at z > 6.5). We find two likely candidates at 
z > 6.5. We attempted Keck/DEIMOS spectroscopy to identify true nature of these candi- 
dates. One candidate has a weak spectral feature which we tentatively identify as Lya at 
z ~ 6.5. Though the colors and photometric redshifts of the second candidate indicate that 
it is likely to be at z > 6.5, the true nature of this candidate is ambiguous. Hence, based 
on our analysis, we predict 1±1 LBG candidates at z > 6.5, brighter than AB~24.5 mag, 
in the GOODS-N field. The expected SFRs of 100-200 M yr" 1 implies that these LBG 
candidates at z > 6.5 are likely to be the progenitors of massive LBGs found at z~5. The 
number density of such luminous LBGs puts strong constraints on the bright-end of the lu- 
minosity function at z > 6.5, consistent within the theoretical limit based on the halo mass 
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function. Deep NIR spectroscopic observations will help to confirm these bright candidates 
as well as shed light on the nature of such luminous LBGs at z > 6.5. 
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Fig. 1. — The Subaru/MOIRCS (dashed) and the CFHT/WIRCam (large box shown with 
dot-dash) field layouts in comparison to the GOODS-N ACS (solid) field. The GOODS-N 
ACS field (a large grid of 40Kx40K pixels) was divided into 17 sections (each 8Kx8K pixels) 
for data handling convenience. The gray solid lines show layout of these sections. 



-26 - 



100 
80 



J 60 
"tn 

| 40 
o 



20 



- 


1 


1 


1 


1 




- 




1 

I J 




H 




Ks 




■ / 




I , .'/ 


I 




i 


1 " 



1 100 



100.0 
0.6 



1253 1427 1625 1850 
Wavelength [nm] 



2107 2400 



Rest-frame Wavelength [nm] 
124.6 155.2 193.3 240.8 



300.0 



0.5 r 



T 



1216 A 
Lyman break 
z~ 



WIRCAM 
Ks-Filter 



WIRCAM 
J-Fllter 



1 J-F?lte,r 




MOIRCS 
H-Filter 



MOIRCS i 
Ks-Filtqr 



- 



997 1241 1547 1927 

Observed Wavelength [nm] 



2400 



Fig. 2.— [Top] The Subaru/MOIRCS J, H, K s (solid) and the CFHT / WIRCam J, K s 
(dashed) filter transmission curves. [Bottom] Location of rest-frame 1216 A Lyman break 
at various redshifts (z > 6.5) and filter throughput curves of z$5o-, J-, H-, i^ s -bands. The 
top axis shows corresponding rest-frame wavelengths at z ~ 7. 
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Fig. 3. — Relative astrometric offsets between the combined -K" s -selected and the GOODS-N 
ACS ^so-selected catalogs. Compact sources with S/N > 20 and FWHM < 1.2" are plotted. 
The histograms show distribution of the offsets in RA and DEC. 
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Fig. 4. — [Top] Comparison between X s -band magnitudes (SExtractor AUTO) obtained 
from the WIRCam and the MOIRCS images. [Bottom] Comparison between K s magnitudes 
obtained from the combined image and K s magnitudes from the 2MASS catalog. We compare 
2MASS magnitudes in a magnitude range between ~14 and ~16 mag (dot-dash vertical lines 
in the bot tom panel) becau se of the non-linearity /selection issues outside this magnitude 



range (see 



Wang etal.ll2010l for details). 
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Fig. 5. — [Left 4-panels] The i^ s -band number counts for each image. The bottom-left panels 
shows MOIRCS f^s-band number counts, while the top-left panel is for the WIRCam. The 
bottom-right panel shows J-band number counts obtained from the combined image, while 
the top-right panel is for i^ s -band. In all panels, the dashed vertical line shows the magnitude 
at which the number counts fall to 50% of their peak value. It is clear that combined images 
are ~0.5 mag deeper than the WIRCam or MOIRCS images. [Right 2-panels] The -fr s -band 
completeness obtained from the simulations. Details of these simulations are discussed in 
text (§ 14. 3p . Top panel shows the completeness from the combined -fT a -band image, while 
the bottom panel is for the combined J. The dashed vertical line shows the magnitude at 
which the recovery rate falls to 50%. The dot-dash curve shows the reliability curve (§ 14.31) 
for both J- and -fT s -selected catalogs. The dot-dash vertical line shows the magnitude at 
which the reliability falls to 50%. 
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Fig. 6. — [Top] Selection of LBGs at z > 6.5 using (z- J) versus (J-K s ) colors . The gray 
shaded region i s the selection region. The red curves are BC03, iKinney et al.l (119961 ) and 



Coleman et al.l ( 119801 ) model colors of low redshift ellipticals, while blue curves (with open 
circles indicating different redshifts starting from the bottom z =6.0, 6.2, 6.4 ...) are BC03 
models of star-forming galaxies with E(B-V)=0, 0.15, 0.30 mag corresponding to three 
differen t curv es. The black filled circles are expected colors of late-type stars based on 



Pickles! ( 119981 ). The black data points (dots) are all objects in the catalog. The filled black 



squares are two zsso-dropout candidates, and the open black squares are 12 other candidates 
selected based on this color criteria before checking their Spitzer magnitudes and photometric 
redshifts. The black data points (dots) in the selected region were excluded by other criteria 
as given in § 15.11 [Bottom] Redshift selection functions at different magnitude bins are shown 
by color-coded histograms. 
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Fig. 7. — Two possible candidates (top panel is candidate A, bottom panel is candidate B) 
at z > 6.5. The coordinates and photometry of these candidates are shown in Table [3j Left 
four stamps are ACS -8435,^06,2775,^850 images at 0.06" pixel -1 , and right four stamps are 
lower resolution J, K s , [3.6], [4.5] images, respectively. The circle shows the position of the 
object and is 1.5" in radius. 
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Fig. 8. — Spectral energy distributions of two z 850 - drop outs. The stellar SED (blue) and a 
galaxy SED (red) are shown for each object. Based on these SEDs, Candidate B is equally 
likely to be a star and a z 85 o-dropout galaxy. If these candidates are z 85 o-dropout galaxies 
then we also show their best-fit stellar ages and stellar masse s based on their SED fits. 
We have used Le PHARE (jArnouts et all 1 19991 : lllbert et all 120061 ) SED/photometric redshift 
code. The TFIT photometry (except SExtractor measured H ma g) is used for SEP fi tting 
because it is much more consistent for multi- wavelength data (e.g.. iDahlen et al.ll2010l ). 
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Fig. 9. — [Top] Keck DEIMOS two dimensional spectrum of candidate B covering 8400- 
10000 A, displayed using histogram equalization, is shown with the object position bracketed 
by horizontal red lines. The gray scale is such that the darker is more positive. [Bottom] 
Corresponding one dimensional spectrum obtain through spec2d reduction software. There 
is a possible weak Lya line detection at ~9090 A in this spectrum, indicated by a dot-dash 
line in one dimensional and marked by a red arrow in the two dimensional spectrum. If 
confirmed, this £ 85 o-dropout candidate could be a LBG at z ~ 6.5. 
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Fig. 10. — Constraints on the bright-end of the luminosity function a t z > 6.5 ( jBouwens et al. 
20081 : besch etal.lhood bold ICapak etallboili bowler etaJboij ) along with our number 
density based on one possible candidate plotted as open red circle. The upper limit is for 
two possible candidates, while lower limit (down arrow) means noz> 6.5 candidates in this 
area at this brightness, and the number density has to be less than this value. The small 
area-deep surveys are not sufficien t to probe the bright-end of the luminosity function. We 
will need wide area surveys (e.g., iTrenti et al.l 1201 ll ; lYan et al.l 1201 ll ; iGrogin et al.l 120111 ) to 
explore these luminous galaxies at these redshifts. 
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Table 1. GOODS-N Near-IR data 



Telescope / Camera 


Filters 


Area Coverage 
(arcmin 2 ) 


Image Pixel Scale 
(arcsec pixel -1 ) 


Exposure 
Time d 


K s Magnitude 
Limit e (mag) 


Typical 
Seeing (") 


Subaru/MOIRCS 


J, H, K s 


109 


0.15° 


-2 hrs 


-24.0 


0.6-0.8 


CFHT /WIRCam 


J, K s 


1040 b 


0.30 


—3 hrs 


-24.0 


0.6-0.8 


Combined 


J, K s 


1040 6 


0.30 


—5 hrs 


-24.5 


0.6-0.8 



a Native pixel scale is 0.117" pixel 1 but to be compatible with other GOODS multi-wavelength data we have 
used 0.15" pixel" 1 . 

b Total WIRCam coverage. For the zsso-dropout selection, we have used only GOODS-N ACS coverage area of 
—169 arcmin 2 . 

c MOIRCS and WIRCam images are combined together to increase the depth in the MOIRCS covered area, which 
is inside the GOODS-N ACS covered region. 

d 80% of the area has at least this exposure time. 

°50% completeness limit (3er) from Figure [5j 
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Table 2. SExtractor Parameters used for /^-selected and J-selected catalogs. The values 
in parenthesis are used for J-selected catalog with other parameters same as _K" s -selected 

catalog. 



Parameter 


Value 


DETECT_MINAREA 


4(2) 


DETECT_THRESH 


2.0 (1.0) 


ANALYSIS_THRESH 


2.0 (1.0) 


FILTER 


Y 


FILTER_NAME 


gauss_3 . _5x5 . conv 


DEBLEND_NTHRESH 


32 


DEBLEND_MINCONT 


0.0001 


CLEAN 


Y 


CLEAN_PARAM 


1.0 (0.3) 


SEEING_FWHM 


0.8 


BACK_SIZE 


80 


BACK_FILTERSIZE 


3 


BACKPHOTCLTYPE 


LOCAL 


BACKPHOTO_THICK 


12 


WEIGHT_TYPE 


MAP_RMS 



Table 3. Photometry of two zsso-dropout candidates. 



ID 


RA 


DEC 


Da 
^435 


v 606 


l 775 


z 850 


J a 


H b 


K s 




m 4.5 




m 8.0 


Stellarity c 


z phot 


A 


189.07061 


62.20895 


>28.0 


31.40 


>28.0 


27.07 


23.22 


23.07 


22.53 


22.18 


22.24 


22.10 


22.04 


0.54 


7.43 










±3.26 




±0.16 


±0.06 


±0.17 


±0.05 


±0.01 


±0.02 


±0.20 


±0.21 






B 


189.08286 


62.15940 


>28.0 


29.90 


28.09 


26.03 


23.94 


24.29 


23.82 


23.62 


24.25 


24.19 


>23.75 


0.13 


6.92 










±0.81 


±0.31 


±0.07 


±0.09 


±0.54 


±0.14 


±0.03 


±0.08 


±0.55 









a Thcse magnitudes are based on TFIT technique |Laidler et al.ll2007t IPahlen et al.lboiol ). 
b These magnitudes are based on SExtractor MAG_AUT0 estimates. 

c Star/Galaxy separation based on SExtractor CLASS_STAR index (l=star) measured in J-band. 



